WKB Analysis of VT- Symmetric Sturm-Liouville problems. II 

Carl M. Bender 1 !!] and Hugh F. Jones^ 

1 Department of Physics, Kings College London, Strand, London WC2R 1LS, UK * 
2 Blackett Laboratory, Imperial College, London SW7 2AZ, UK 

In a previous paper it was shown that a one-turning-point WKB approximation 
gives an accurate picture of the spectrum of certain non-Hermitian PT-symmetric 
Hamiltonians on a finite interval with Dirichlet boundary conditions. Botentials to 
which this analysis applies include the linear potential V = igx and the sinusoidal 
potential V = igsin(ax). However, the one-turning-point analysis fails to give the 
full structure of the spectrum for the cubic potential V = igx 5 , and in particular it 
£X fails to reproduce the critical points at which two real eigenvalues merge and become 

a complex-conjugate pair. The present paper extends the method to cases where the 
WKB path goes through a pair of turning points. The extended method gives an 
extremely accurate approximation to the spectrum of V = igx 3 , and more generally 
it works for potentials of the form V = igx : . When applied to potentials with 
half-integral powers of x, the method again works well for one sign of the coupling, 
namely that for which the turning points lie on the first sheet in the lower-half plane. 
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The study of PT-symmetric Hamiltonians [l|-|6| has largely concentrated on eigenvalue 
problems in which the wave function is required to vanish at infinity in various Stokes wedges. 
However, some work has also been done on Sturm-Liouville problems defined on a finite 
interval (7| , where the wave- function is required to vanish at the end points. Such problems, 
especially with the potential V(x) = igx, are relevant to a number of physical situations [8j- 
[10|. In a previous paper [11| we showed that for the imaginary linear potential and for 
potentials of the form V(x) = igsin(ax), the WKB approximation with the WKB path 
passing through a single turning point gives a spectacularly good picture of the spectrum, 
even for low energies. In particular, it correctly locates the critical points, where pairs of real 
eigenvalues merge and become complex conjugates of one other. However, for the potential 
V(x) = igx 3 the method failed to produce critical points, and the question of whether a 
different WKB path could do so was not addressed. For the igx 3 potential there are three 
turning points, and an alternative path passes through a pair of PT-symmetric turning 
points. 

In the present paper we develop the formalism for such a path and show in Sec. II that 
the two-turning-point WKB approximation gives a surprisingly simple secular equation. In 
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Sec. Ill we demonstrate that this approximation accurately reproduces the full structure of 
the spectrum for V(x) = igx 3 . The method is also shown to give good results for higher- 
power potentials of the form V(x) = igx 2N+1 , where N is integral. When the method is 
applied to potentials with half-integral powers of x, the double- valued nature of the potential 
must be taken into account. We reproduce the spectrum when there exists a simple path 
that passes through the turning points and remains on the first Riemann sheet. In Sec. IV 
we give some brief concluding remarks. 



II. WKB CALCULATION OF EIGENVALUES AND CRITICAL POINTS 

We consider the eigenvalue equation 

- ij)"{x) + V(x)if>(x) = \ip(x) (1) 

defined on the interval [—L,L], with Dirichlet boundary conditions ip(±L) = 0. In what 
follows we take L = 1 without loss of generality because L = 1 can be obtained by rescal- 
ing for the monomial potentials V(x) = —g{ix) M that we consider. The equation and the 
boundary conditions are PT symmetric, so we expect real eigenvalues when the PT sym- 
metry is unbroken; that is, when the wave functions are themselves PT symmetric. When 
the PT symmetry is broken, the eigenvalues occur in complex-conjugate pairs. 
With \ = ag the equation becomes 

ij"(x)+gQ(x)ij(x)=0, (2) 

where Q(x) = a — V(x) and a is considered as a fixed parameter. For large g the WKB 
approximation to ip(x) has the form 



ijj(x) 



(gQ) 1/4 



A exp (i ds \f 'gQ(s) J + B exp I — i / ds \J gQ{s) J 



(3) 



However, the path connecting —1 to +1 is not specified. If it goes directly from —1 to +1 
without passing through a turning point, ip(x) has the same form along the path. Then, by 
imposing the boundary conditions ^(±1) = 0, one obtains the standard no-turning-point 
approximation 

sin I T = 0, (4) 

where It = J_ 1 ds \/gQ{s). The eigenvalues are then obtained from this secular equation, 
which always gives real eigenvalues and predicts no degeneracies (critical points) in the 
spectrum. 

When, as in Ref . [ri| , the path goes through a single turning point on the imaginary axis 
at xq = ib, the wave function has different coefficients to the left and to the right of the 
turning point. These coefficients are determined by matching asymptotically the separate 
WKB approximations to the Airy functions that approximate the solution to the Schrodinger 
equation in the vicinity of the turning point. In that case the resulting secular equation is 
found to be 

sinJ T + |e A = 0, (5) 

where A = 2Im/ + with I + = j ds \/ gQ(s). Because of PT symmetry, the total integral 
It can be written as It = 2Re/ + . 



It is the presence of the second term in this new secular equation that gives the interesting 
structure and critical points of the spectrum. When A is real, A is real and intrinsically large 
in magnitude. When A is negative, the second term is generically very small, so that the 
equation essentially reduces to (@J; when A is positive, the second term is generically large 
and cannot be balanced by the first term unless A becomes complex. The critical points 
where this change of behaviour occurs are given by A <^C 1. 

In the case of the potential V(x) = igx 3 there are three turning points at x = ib(l, u>, co 2 ), 
where b = a 1//3 and u> = e 2l7r//3 . Numerical calculations [ll| show a series of critical points, 
which the one-turning-point WKB approximation fails to reproduce. We are therefore mo- 
tivated to develop the formalism for a two-turning-point approximation, where the path of 
integration passes through the turning points xl = ibu and Xr = ibu 2 . In this case we need 
three different WKB approximations to the wave function: ipi{x) for % between —1 and Xl, 
iPm(x) f° r x between xl and Xr, and iJJr(x) for x between xr and 1. These approximations 
take the form 



ip L (x) 
iPm(x) 

i>B.{x) 



(gQ) 1/4 
l 

WW 4 
l 

WW 1 



L\ exp I i / ds \/gQ{s) j + L 2 exp ( — i / ds \fgQ{s 
M\ exp I i / ds \J gQ(s) J + M 2 exp \—i ds \f 'gQ(s) 
Ri exp I i / ds \JgQ{s) J + R 2 exp \—i ds \/gQ(s) 



(6) 



We will not need tPr(x) explicitly, as our tactic will be to work with i/)l(x) and iPm{x) in 

the left half-plane and to enforce PT symmetry of the wave function on the imaginary axis. 

The first constraint on the coefficients is that iPl{%) should vanish at x — — 1, which gives 

the equation 



L 2 



-2ilr 



(7) 



where II = /_f dsy/gQ(s). The second constraint comes from requiring that the wave 
function be PT symmetric, namely, that ip'(x)/ij)(x) be pure imaginary at an arbitrary 
point x = —ia on the imaginary axis. After some algebra this leads to the equation 



MlMl - e -^M 



M 2 M{ 



(8) 



where I M = f** ds^/gQ(s). 

Finally, we match the two WKB approximations to the superposition of Airy functions 



ij) A (x) = KiAi(y) + K 2 M{u 2 y), 



(9) 



with y = (x — xl)/c, that approximates the solution to the Schrodinger equation in the 
vicinity of xl- The coefficient c is given by c = , ye~ ie ^ 3 , where 7 = [b/(3ga)] and 6 = 5ti/6. 
Because g is large, c is small. Hence, y is large, and the Airy functions have two distinct 
asymptotic behaviours depending on the argument of y, namely 



Ai(y) 



1 



2 v ^F 2 / 1 /4 



-2,3/2 



(|arg y\ < n) 



Ai(y) 



2v/i?/ 1/4 



-2„3/2 



e 3y + iez y 



2«3/ 2 



|arg j/ = vrj 



(10) 



In the neighbourhood of xl the two WKB approximations can be further approximated as 



iMs) oc -L (M ie l^ 3/2 + M 2 e-i^ 3/2 



(11) 



In the complex-?/ plane the matches must be made along the Stokes lines where the wave 
function is purely oscillatory. These Stokes lines are indicated in the left panel of Fig. [TJ The 
rotation angle between x and y is 9/3, which is approximately 7r/3, and the paths have to 
be adjusted so that y and u 2 y lie on Stokes lines. The required directions in the complex- £ 
plane are indicated in the right panel of Fig. [TJ For the matching to i/)m(%), both y and 





FIG. 1: Matching paths in the complex y-plane (left panel) and x-plane (right panel). 



u 2 y lie away from the negative-?/ axis. Thus, in this case we can use the first asymptotic 
approximation of (TTOj) . This results in the matching 

M ie l yi/2 + M 2 e-^ /2 oc K ie -^ /2 + K,< ' ; ( '< $" 
and leads to the condition 



|„3/2 



(12) 



M 1 K 2 



in/6 



(13) 



M 2 K t 

For the matching to iPl(x), u 2 y again lies away from the negative-;?/ axis, but y lies along 
that axis. In that case we use the second approximation of (fTUj) for Ai(y). The resulting 
match is 

(14) 



-r -lifi/ 2 -r lifl/ 2 -r^ / -2„3/2 . 2 3/2 \ T ^ ,• / fi 2 3/2 

Lie 'i y + L 2 ez y oc K x [e z y +ie* y + K 2 e w b e^ y , 



which gives the condition 



^ - 7 + ^ P «/6 



(15) 



Combining flTJ, ©, (USD and ([15]), we obtain 

7 J- p2iI L 

e~ 2ilM = !-«/> » ( 16 ) 

— « + e £ 

which results in the surprisingly simple secular equation 

sin It + e cos Im = 0, (17) 

where A = 21iu.Il. 

The secular equation ( TTTj) is our main result. It was derived explicitly for the case V(x) = 
igx 3 , but in fact it applies without change for potentials of the form V(x) = —g(ix) 2N+1 
when the two relevant turning points lie in the lower-half plane. If the relevant turning 
points lie in the upper-half plane, the sign of A in ( 117)) should be reversed. When N is an 
integer, the spectrum is the same for g negative because changing the sign of g is equivalent 
to changing the sign of i. Otherwise, the spectrum is asymmetric in g. 

A few words about the role of A are in order. As already mentioned, it is the change 
in sign of A that produces the interesting structure in the spectrum and determines the 
location of the critical points. The significance of A = is that when A vanishes, II is 
purely real and there is a purely oscillatory path for the WKB wave function in going from 
x = — 1 to x = xl- There is a path connecting xl and xr along which the wave function 
is purely oscillatory, so the wave function is oscillatory along the entire path from -1 to 
1. If there is no oscillatory path starting from x — — 1 that passes through xl, the WKB 
approximation fails to produce critical points. These considerations guide us in the choice 
of the appropriate turning points. 

III. NUMERICAL CALCULATIONS 

For V = igx 3 the two relevant turning points are shown circled in Fig. [2] for a = 0.5. 
We also show the Stokes lines emanating from xl and xr along which the exponent in the 
WKB approximation is purely imaginary. There is a purely oscillatory path between xl and 
Xr, and the continuations towards ±1 almost pass through those points. They do so when 
A = 0. 

Figure [3] shows the result of (ITTj) for the real spectrum of the potential V = igx 3 along 
with the numerical results obtained by a shooting method. As can be seen, the WKB 
approximation reproduces the numerical result very accurately. The blue dots extending a 
little way beyond a critical point represent the common real part of the complex conjugate 
pair of eigenvalues. The spectrum for g negative is not shown because it is identical to that 
for g positive. 

Figures H] and show the corresponding results for the potential V = —igx 5 . Because 
the relevant turning points are above the real-x axis, we must now use e _A in (ITT)) . Again, 
the WKB approximation reproduces the numerical result extremely accurately. A general 
feature of the spectrum is that as the power M of x increases, the interesting structure 
occurs for higher values of g. This is readily understood because the problem is initially 
posed in the range — 1 < x < 1. 

It seems clear that this WKB approximation will work equally well for potentials of the 
form V = —(igx) 2N+l , where N is integral. The turning points chosen are probably those 
closest to the x-axis, but in any case they must be a pair for which the Stokes lines from xl 
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FIG. 2: Turning points for the potential V = igx 3 with a = 0.5. The two turning points used in 
the WKB approximation are circled and marked xl and xr. Also shown (dashed green lines) are 
the Stokes lines emanating from xl and xr. 
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FIG. 3: WKB approximation (red line) for the real energy levels of the potential V(x) = igx 3 
compared with the numerical results (blue dots). The WKB approximation and the numerical 
results are virtually indistinguishable. 



and xr pass through x = ±1 for some value of a. This corresponds to A = and the onset 
of structure in the spectrum, resulting in critical points. 

We have also examined the effectiveness of the WKB approximation for half-integral PT- 
symmetric potentials. Specifically we have looked at the cases V = —g(ix) M for M = 1/2, 
3/2 and 5/2. Now we have to differentiate between the cases g positive and g negative. 
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FIG. 4: Turning points for the potential V = —igx 5 with a = 0.5. The two turning points used in 
the WKB approximation are circled and marked xj_, and xr. Also shown (dashed green lines) are 
the Stokes lines emanating from xl and xr. 
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FIG. 5: WKB approximation (red line) for the real energy levels of the potential V(x) = —igx 5 
compared with the numerical results (blue dots). As in Fig. [3l the WKB approximation and the 
numerical results are virtually indistinguishable. 



Starting with V = — g(ix) 1//2 , the potential is a two-sheeted function with a cut that 
we take to lie along the positive imaginary axis. There is a single turning point, which is 
located at x — —ia 2 , but it is on different Riemann sheets depending on the sign of g. For 
g < the turning point is on the first sheet; however, for g > it is on the second sheet. 
It is natural to take the end points x = ±1 to be on the first Riemann sheet. Then, for 



g < there is a straightforward WKB path through the turning point. On the other hand, 
for g > the path would have to pass through the cut twice in order to pass through the 
turning point. We do not know how to implement such a path, so for g > we simply use 
the straightforward no-turning-point formula of (jl]). By its nature this latter approximation 
is rather smooth and does not reproduce the critical points that occur in the numerical 
solution. The two WKB approximations and the numerical data are presented in Fig. |6j 
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FIG. 6: WKB approximations (red line) for the real energy levels of the potential V(x) = —g(ix)2 
compared with the numerical results (blue dots). For g > we use (jl]), while for g < we use ([5]) 
with e A replaced by e~ A . See text. 



Going on to M = 3/2, the turning points are shown in Fig. [7J There is a single turning 
point, marked Xo, at x — —ia 2 ^ 3 , and a pair of turning points, xl and xr, in the upper half 
plane. For g < the single turning point Xq is on the first sheet, and the appropriate path 
is one going through that point, giving the WKB approximation of (J3J). As can be seen 
from Fig. 8, the approximation is extremely good. For g > 0, the single turning point is on 
the second sheet, while xl and xr are on the first sheet. However, there is no Stokes lines 
between them because of the presence of the cut. This explains why the WKB approximation 
of (fT7|) with e A replaced by e~ A fails to reproduce completely the numerical structure of 
the spectrum. Finally we turn to M — 5/2. The relevant turning points are shown in 
Fig. [9j For g > the turning points xl and xr are on the first sheet, and the appropriate 
path is one going through them, giving the extremely accurate WKB approximation of (|17|) . 
For g < 0, the upper pair of turning points are on the first sheet, but again there is no 
Stokes lines between them because of the presence of the cut and consequently the WKB 
approximation fails to reproduce completely the numerical structure of the spectrum. 

IV. CONCLUDING REMARKS 



We have derived a very simple secular equation for the real eigenvalues of PT-symmetric 
Sturm-Liouville problems where the potential is of the form V = igx 2N+1 with N integral. 
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FIG. 7: Turning points for the potential V(x) = —g(ix) 3 ' 2 , for a = 0.7. 
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FIG. 8: WKB approximations (red line) for the real energy levels of the potential V(x) = —g(ix) 3 ' 2 
compared with the numerical results (blue dots) . For g < we use ([5]) , while for g > we use (fTTj) 
with e replaced by e . This fails to reproduce the spectrum because of complications due to 
the cut. 



This equation is derived using the WKB approximation for a path that passes through the 
pair of turning points lying nearest to the x-axis, and supplements the equation previously 
derived for a path passing through a single turning point. For the above potentials the 
spectrum is symmetric in g, so that it is only necessary to consider one sign of g. The 
spectra so obtained are extremely accurate and show all the interesting structure, including 
the critical points where two eigenvalues merge and become complex conjugates. 

We also considered the extension of the secular equation to potentials with half-integral 
powers. In those cases the situation is complicated by the presence of a cut, and a straight- 
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FIG. 9: Turning points for the potential V(x) = —g(ix) 5 ' 2 , for a = 0.7. 
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FIG. 10: WKB approximations (red line) for the real energy levels of the potential V{x) = 
—g(ix) 3 ' 2 compared with the numerical results (blue dots). For g > we use (|17p with the 
turning points xl and xr, while for g < we use the upper pair of turning points. This latter 
approximation fails to reproduce the spectrum because of complications due to the cut. 



forward application is only possible for the particular sign of g where the relevant critical 
point or points lie in the lower-half plane. Here the results are again extremely accurate. 
For the other sign of g the method fails to reproduce the complete structure of the spectrum, 
and an extension of the method is needed to deal with the situation when either the turning 
points are on a different sheet from the end-points ±1 or the path between them is impeded 
by the cut. 
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As mentioned in the Introduction, this problem with the linear potential V = igx occurs 
in a number of physical situations, and one may speculate that the higher-power potentials 
considered in this paper might also have some physical relevance. The most likely application 
is in the context of hydrodynamics, discussed in Ref. [10|, where the potential V(x) is 
associated with the velocity profile q(x) of the fluid. 
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